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Abstract 

We report on a new class of electromagnetically-driven fluid interface instability. Using the 
optical radiation pressure of a cw laser to bend a very soft near-critical liquid-liquid interface, we 
show that it becomes unstable for sufficiently large beam power P, leading to the formation of 
a stationary beam-centered liquid micro-jet. We explore the behavior of the instability onset by 
tuning the interface softness with temperature and varying the size of the exciting beam. The 
instability mechanism is experimentally demonstrated. It simply relies on total reflection of light 
at the deformed interface whose condition provides the universal scaling relation for the onset Ps 
of the instability. 

PACS numbers: 47.20.Ma, 42.50.Vk, 82.70.Kj, 42.25. Gy, 47.27.Wg 



1 



Deformations and break-up of fluid interfaces under an applied field play a significant 



r „,e in scientific and technic endeavor, Since the ear* wo rks of Ze.env Q and fn rt ne r 
investigations of Taylor ^, interface instabilities driven by electric fields are the most famil- 
iar examples. They represent nowadays the corner stone of numerous industrial processes 
as different as electro-spraying Q], electrospinning of polymer fibers Q or surface relief pat- 
terning From the fundamental point of view they also illustrate a simple and fascinating 
example of behavior leading to finite time singularity gj]. A uniform magnetic field can as 
well destabilize fluid interfaces and form well-organized peak structures |2| or elongate mag- 
netic droplets though no liquid ejection occurs in that case. These deformations, as well 



as those driven by the acoustic radiation pressure [9], were essentially used to characterize 
in a non contact manner, the mechanical properties of fluid interfaces |~ul [~Lll|. 

Recently, with the increasing importance of optical forces in emerging nano/bio- 
technologies, this strategy has been successfully extended to the optical deformation of 
soft materials or micro-obi ets and brought new insights into their viscoelastic properties 

nnn 

|12L 111 LLJ]. However, as the bending of fluid surfaces by the optical radiation pressure is 
usually weak, further developments require larger beam intensities with smaller disturbing 
effects. Interface disruption under the applied field could constitute a major drawback for 
this promising method and should therefore be investigated. Yet, even if non linear elonga- 
tions of red blood cells Q| and droplets disruption [3] were already shown, to the best of 
our knowledge interface instability driven by the optical radiation pressure has never been 
observed nor discussed in detail. 

Using very soft transparent liquid interfaces to enhance optical radiation pressure effects, 
we show in the present Letter that, contrary to expectations, the hump induced by a con- 
tinuous laser wave becomes unstable for sufficiently large beam powers. To analyze this new 
opto-hydrodynamic instability, we explore the universal behavior of its onset on near-critical 
liquid-liquid interfaces by varying the interface softness with a temperature scanning, and 
by changing the size of the exciting beam. When properly rescaled, the dispersion in mea- 
sured onsets reduces to a single master behavior which is retrieved theoretically from simple 
arguments. In fact we demonstrate that the total reflection condition of light at the inter- 
face defines an onset in power which agrees with the observed scaling law. Finally we show 
that the optically-driven interface instability leads to the formation of a stationary beam- 
centered liquid micro-jet emitting droplets, which anticipates the bases for new applications 



2 



in microfluidics and liquid micro-spraying. 

Experiments. — Opto-hydrodynamic instabilities of fluid interfaces are realized in a qua- 
ternary liquid mixture made of toluene, sodium dodecyl sulfate (SDS), n-butanol and water. 
At room temperature we prepare a water-in-oil micellar phase of microemulsion whose com- 
position in weight fraction is 70 % toluene, 4 % SDS, 17% n-butanol and 9% water. Close 
to Tc, where Tc = 35°C is the liquid-liquid critical temperature, the mixture belongs to the 
universality class (d=3,n=l) of the Ising model For a temperature T > Tc it separates 
in two micellar phases of different concentrations $i and $2- The use of these mixtures was 
motivated by the fact that significant interface deformations by optical radiation, without 
nonlinear propagation effect or disturbing thermal coupling require weak surface tension 
and buoyancy. Indeed both effects vanishes when the critical point is neared respectively 
with the critical exponents 2u = 1.26 and (3 = 0.325. Then the surface tension a of phase- 
separated supramolecular liquids is intrinsically small. For example from a = o"o(^ — l) 2i/ , 
where <7o = 1.04 10~ A N/m in our case, we get a ~ 10~ 7 N/m at (T — Tc) = 1.5K a value 
typically 10 6 times smaller than that of the water free surface at room temperature. 

The experimental configuration is schematically presented in Fig. la. The mixture 
is enclosed in a fused-quartz cell of path length 1=2 mm that is thermally-controlled at 
the temperature T above T c (temperature accuracy: AT ^ 0.05K). Since density (resp. 
index of refraction) of water is larger (resp. smaller) than that of toluene, the micellar 
phase of larger concentration $i is located below the low micellar concentration phase $2, 
while its refractive index n\ is smaller than 72,2 of $2- The difference in density is given by 
P\ — P2 = (Ap)o(^ — I) 13 with (Ap) ~ 285 kg.m~ 3 , and the Clausius-Mossotti relation close 
to the critical point leads to n 2 — n\ = (§^)r(P2 - Pi) with (| 2 )r ~ -1.22 10~ 4 m 3 .kg~ 1 . 
The bending of the liquid-liquid meniscus is driven by a linearly polarized TEM 00 cw Ar + 
laser (wavelength in vacuum, A = 5145 A) propagating vertically downward from $2 to 
$1. The beam is focused on the interface by a 10X microscope lens to ensure a cylindrical 
symmetry of the intensity profile around the meniscus. We adjust also the beam-waist ujq 
by changing the distance between a first lens (focal length f=l m) and the 10X objective. 

As already demonstrated, the optical radiation pressure originates from light momentum 
discontinuity at the interface. It always pushes normally to this interface the fluid with 
larger refractive index into that with smaller refractive index, regardless of the direction of 
propagation [r^ . As a result, radiation pressure acts downwards in our case (see Fig. 1), and 
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the height h(r) of the induced deformation is determined by a balance with the buoyancy 
and Laplace restoring forces. Then, at steady state h(r) is described by ["b^ : 

/ n , / \ Id/ rh (r) \ „, s 

to _ ft)9Mr) _ ff __(_U_) = n W (i) 

Primes denote derivatives versus r and g is the gravity constant. In the general case, the 

optical radiation pressure Il(r) is given by (lf^ : 

n 2 cos 2 6 2 \ tan 2 ] . . , 

n(r = 1 + — —TI{r) 2 

c L tan d\ J 

where J(r) = ^ ex P(-^r-) is the intensity of the incident wave, P the beam power and 
c is the light velocity in vacuum. R and T are the Fresnel coefficients of reflection and 
transmission in energy. 9 2 and 6\ are respectively the incident and the refraction angles 
versus the interface normal and are given by cos #2 = ; 1 , „ : and the Snell relation. 

/ y/l+h' (r) 

Results. — For small bending curvatures, i.e. for ft (r) < 1, deformations are function of 
an optical Bond number Bo defined as the ratio of buoyancy over the Laplace force 



Bo = (j£) 2 , where lc = y ( Pl - P2 ) g 18 the capillary length. In the present analysis, we used 
values of ujq and (T — Tc) ranging respectively from 3.5 to 11 fim and 1.5 to 15K in order to 
work in the Bo < 1 regime (we have experimentally 0.006 < Bo < 0.54) and reproduce the 
case of classical liquid free surfaces. A typical evolution of the induced deformations versus 
incident beam power P is illustrated in Fig. la-lh for T — T c = 15K and u = 3.5/im, i.e. 
for the smallest Bo investigated. Fig. 2 depicts also the variation of h(r = 0) versus P for 
the two extreme values of our Bo number range; experimental errors are ^ ~ < 5% . 
As expected, h(r = 0) is linear in P as far as h (r) <^ 1, i.e. typically for P < 970 mW in 



the example presented (Fig. la-le). T 
with the linear behavior predicted 



le measured slope is moreover in complete accordance 
[. Then, with further increase in P, h(r = 0) rapidly 
deviates from the linear regime (970 < P < 1170 mW) and suddenly diverges when P 
reaches a power onset P$ = 1230 mW (see the transition from Fig. lg to Fig. lh). Fig. 
2 shows that this behavior is observed over a large range of the parameter space and that 
instability occurs even more rapidly (Ps = 370 mW) when the critical point is neared 
(T — Tc = 1.5K ), despite an increase in beam- waist (u = 11 fim). 

Consequently, since Ps vary noticeably in the range of Bo investigated, we analyzed its 
behavior over large variations in ojq and (T — Tc). At first, as illustrated in Fig. 3 over one 
order of magnitude in (T — Tc), Ps varies linearly versus u : the instability onset decreases 



with focusing. Moreover results show that the slope of Ps versus too also decreases when the 
critical point is neared. This behavior is much better illustrated in the Inset of Fig. 3 for 



vl.01±0.05 



the smallest beam-waist used loq = 3.5/xm. A power law fit leads to Ps oc (T — Tc j 

Discussion. — While interface deformations do not depend on the direction of propagation 
at low laser excitations, this is no more true at higher field strengths. Instead of a liquid 
jet, huge tethered deformations are induced by an up ward propagating beam, but they 
remain always stable, despite surprising aspect ratios |2£j| . A simple physical mechanism that 
differentiates both cases is the possibility to reach total-internal reflection when the beam 
propagates from $ 2 to $1, i.e. from the large to the low refractive index phase. Furthermore, 
as it can be seen from Eq. (2), the transfer of momentum at the interface is maximum under 
total reflection condition because R = 1 and T = in that case. A comparison with the 
momentum transfer at normal incidence shows that the radiation pressure is thus enhanced 
by a factor (1 + ~ 4 (ni ~ n 2 in our case). This factor is obviously underestimated 
because secondary reflections and even self-guiding of light inside the deformation will also 
contribute to a further increase. Then, let us express the total reflection condition: 

h'(r) n-i , , 

sinfl 2 = -=y > — 3 

^/l + h'{r) 2 n 2 v ; 

Since experiments are realized in the Bo < 1 regime, buoyancy does not play any important 
role and can thus be neglected. On the other hand, we expand the expression of the radiation 
pressure given by Eq. (2) to first order in n 2 — ri\ because $ 2 and $i are coexisting phases 
of close composition due to the vicinity of the critical point. Eq. (1) becomes: 

r dr \ y/i + /i'( r ) 2 / c Vnj + ni/Trwg V uj$ J 
Integration of Eq. (4) gives directly the expression of sin 9 2 and Eq. (3) leads to: 



n 2 fn 2 — n\\ P 1 



-2r 



2 



exp( 



OJ, 



(5) 

n 2 



c \n 2 + n\J ira r 

Total reflection will then be locally reached for = 1.121, which renders the left hand 
term of Eq. (5) maximum. This condition defines an onset in power 

„ 1.121 7r ni A, n>r\ oc 

P> p — 1 + — uj 6 

0.715 V?n 2 V n 2 J n 2 - m 

which is linear in u and behaves as (T — T c ) 2v ~' 13 , i.e. with the exponent 2v — (3 = 0.935, 
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in very close agreement with the critical behavior of Ps- In looking for scaling, Eq. (6) 
also provides a natural way for presenting the instability onsets. This is illustrated on 
Fig. 4 where the whole set of data is brought onto a single master curve. As expected, 
a linear behavior of Ps versus is observed and a linear fit leads to PsiW) = (6.3 ± 

0.3)10~ 6 [ ra2 °^ n J Uo(fim) while the predicted slope is 6.9 1CT 6 . Optohydrodynamic interface 
instability is then clearly triggered by the total reflection of the exciting beam inside the 
induced deformation. Hence, the brightness of the liquid jets (see Fig. li,j) confirms that 
they behave as liquid waveguides. 

Finally, since matter with large refractive index is always trapped in high beam intensity 
regions by optical forces j^jj, the optical destabilization of liquid interfaces leads to the 
formation of a stationary beam-centered jet which eventually breaks into a spray of nearly 
monodisperse droplets (Fig. li). Bearing in mind that our system has ultralow surface 
tension and that the density difference between the phases is also weak, it can be calculated 
that the radius of the induced jet is much smaller than the typical size for which fluctuation 
effects become important for pinching [3,0, Q|- The break-up mechanism, which deserves 
further investigation, is very robust and convenient because inclined jets can also be easily 
created by tilting the laser beam, as shown on Fig. lj. 

In conclusion, we have presented a new electromagnetic instability mechanism of fluid 
interfaces driven by the optical radiation pressure of a cw laser wave. It opens the way 
to optical streaming of fluid interfaces, even if laser pulses could be necessary for stiffer 
interfaces The stability of the resulting jet, as well as the regularity of the produced 
micro-droplets, offers also new perspectives for optical guiding of microfluidic flows. 
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FIG. 1: (a-h) Variation of the optical interface bending versus beam power until instability when 
(T - T c ) = 15K and u = : (a) P = 310 mW, (b) P = 500 mW, (c) P = 620 mW, (d) P = 

740 mW, (e) P = 920 mW, (f) P = 1110 mW, (g) P = 1170 mW , and (h) P = P s = 1230 mW . For 
P = Ps the interface becomes unstable leading to the formation of a stationary liquid jet similar 
to that illustrated in the overview (i) for (T — Tc) = 6K, ujq = 3.5fj,m and P = Ps = 490 mW. (j) 
Inclined liquid jet obtained for the same experimental conditions as (i) but for P = 770 mW > Ps- 
The laser beam is represented by the arrows, and the height of (i) and (j) is 1mm. 
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FIG. 2: Evolution of the height of the deformation h(r = 0) versus P for the extreme values of the 
Bond number range investigated. Parameters are respectively: (•) T — Tc = 15K and ujq = 3.5/xm 
corresponding to the profiles represented in Figs, (la-lh) and (■) T — Tc = 1.5K and ujq = ll/^m. 
Full lines are linear fits of the weak deformation regime. Broken lines indicate the onset Ps above 
which the interface is unstable. 
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FIG. 3: Variation of the instability onset Ps versus beam waist lvq as a function of T — Tc. 
Parameters are: (T - Tc) = 1.5K(»), SK(U), 5K(+), 6K (o), 8K(A), 12K(ffl), and 15K(f). Lines 
represent linear fits. Inset: Log-log plot of the variation of Ps versus T — Tc for ojq = 3.5/xm. The 
full line is a power-law fit giving P s oc (T - r c )i-Oi±o.05_ 
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FIG. 4: Rescaling of the set of dispersed data presented in Fig. 3 according to the behavior predicted 
by Eq. (6). Symbols are the same. The full line represents a linear fit. 
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